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Abstract. We present the results of Monte Carlo simulation studies of the capability of the 
proposed NEMO km 3 telescope to detect TeV muon neutrinos from Galactic microquasars. 
In particular we determined the detector sensitivity to each known microquasar, optimizing 
the event selection in order to reject the atmospheric background. We also determined the 
expected number of source and background events surviving the selection. 
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1. Introduction 

Microquasars are Galactic X-ray binary 
systems which exhibit relati vistic jets, ob - 
served in the radio band (IChatvl I2005I) . 
Several authors propose microquasar jets as 
sites of acceleration of charged particles up 
to energies of about 10 16 eV, and of high 
energy neutrino production. According to 
present models, neutr inos could be pro- 
duced both in py dLevinson & Waxmanl 
1200 It iDistefano etaTI 120021) and p p in- 
teraction scenarions dBednarekl 120051; 



Aharo nian et al. I 120051; ICh ristiansen et al. 
20051: iRomero & Orellanall2005l) . 

The aim of this paper is to study the possi- 
bility to detect neutrinos from known micro- 
quasar s with the proposed NEMO-km 3 tele- 
scope (Migneco et al 1120061) . In particular, for 
each microquasar we calculated the expected 
sensitivity, optimizing the event selection to re- 
ject the atmospheric neutrino and muon back- 
grounds. 
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We also calculated, according to the 
present theoretical models, the expected num- 
ber of microquasar events that survive the 
event selection. We compared the expected 
source signal with the remaining background 
to establish its detectability. 

The NEMO proposal is made in the con- 
text of th KM3 prpj ect for a km 3 detector in t he 
Mediterranean sea dKM3NeT website 1120071) . 

2. The NEMO km 3 detector 

The NEMO-km 3 telescope, simulated in this 
work, is a square array of 9 x 9 towers with a 
distance between towers of 140 m. In this con- 
figuration each tower hosts 72 PMTs (with a 
diameter of 10"), namely 5832 PMTs for the 
whole detector with a total geometrical vol- 
ume of ~ 0.9 km 3 . We considered an 18 storey 
tower; each storey is made of a 20 m long 
beam structure hosting two optical modules 
(one downlooking and one looking horizon- 
tally) at each end (4 OMs per storey). The ver- 
tical distance between storeys is 40 m. A spac- 
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ing of 150 m is added at the base of the tower, 
between the anchor and the lowermost storey. 

The detector response is simulated us- 
ing the simulation codes dev eloped by the 
ANTARES Collabor ation dAmram et al. I 
l2002t lBecherinill2006l) . modi fied for a km 3 
telescope dAiello et all 120071) . In the sim- 
ulation codes, the light absorption length, 
measured in the site of Capo Passe ro (L a a 68 
m at 440 nm dRiccobene 1 120071) ). is taken 
into account. Once the sample of PMT hits 
is generated, spurious PMT hits, due to the 
underwater optical noise ( 4() K decay), are 
introduced, with a rate of 30 kHz for 10" 
PMTs, corresponding to the average value 
measured in the Capo Passero site. 



to the decrease of the time per day spent by the 
source below the Astronomical Horizon (with 
respect to the latitude of the Capo Passero site). 
In Tab. Q] we present only results concerning 
microquasars which may be observed by a tele- 
scope located in the Capo Passero site (i.e. with 
declination 5^ + 54°). 



Table 1. Detector sensitivity to neutrinos from 
microquasars: the sensitivity / Vi 9o (expressed 
in erg/cm 2 s) is calculated for an s~ 2 neutrino 
spectrum in the energy range 1-100 TeV, for a 
detector livetime of 1 year. The corresponding 
values of the search bin angular radius r\, m and 
the source declination 6 are also given, both ex- 
pressed in degrees. 



3. Calculation of the sensitivity 

The detector sensitivity was calculated ac- 
cording to the [feldman & Cousins! (1 19981) ap- 
proach. The 90% c.l. sensitivity to a neutrino 
flux coming from a microquasar is given by 



_ /%(j0 f th 
/v,9() - —j^rJy . 



(1) 



where ju 90 (b) is the 90% c.l. average upper limit 
for an expected background (atmospheric neu- 
trinos + muons) with known mean value b and 
no true signal dFeldman & Cousinsll 19981) . f' v h 
is the theoretical neutrino energy flux from a 
given microquasar that induces a mean sig- 
nal TV™. During the calculation, an event se- 
lection is applied in orde r to optimize t he sen - 
sitivity, as described in Aiel lo et aTl d2007l) . 
Detailed calculation of the sensitivity for the 
proposed NEMO km 3 telescope to a generic 
p oint-like muon neu trino source are presented 
in lDistefano I (120071) . 

The detector sensitivity was calculated for 
a livetime of 1 year, simulating a neutrino flux 
with spectral index T = 2 in the energy range 1 
- 100 TeV. The study was carried out for each 
microquasar, since the sensitivity is a function 
of the source astronomical declination. Results 
are given in Tab.[T]/ v ,9o- Fig.Q]shows the detec- 
tor sensitivity for the studied microquasars as a 
function of the declination; the sensitivity flux 
limit increases with increasing declination, due 



Source name 


run 


./ 


v.90 




6 


Steady Sources 


LS 5039 


0.9 


6.5 


io- 


11 


-14.85 


Scorpius X- 1 


0.7 


5.8 


io- 


1 1 


-15.64 


SS433 


0.8 


5.7 


io- 


1 ] 


4.98 


GX 339-4 


0.5 


4.7 


io- 


1 1 


-48.79 


Cygnus X-l 


0.7 


9.0 


io- 


1 ] 


35.20 


Bursting Sources 


XTE J1748-288 


0.9 


5.4 


io- 


11 


-28.47 


Cygnus X-3 


0.8 


1.1 


io- 


10 


40.95 


GRO J1655-40 


0.7 


5.2 


io- 


1 1 


-39.85 


GRS 1915+105 


0.8 


7.4 


io- 


1 1 


10.86 


Circinus X- 1 


0.9 


4.2 


io- 


1 ] 


-56.99 


XTE J 1550-564 


0.9 


4.4 


io- 


1 1 


-56.48 


V4641 Sgr 


0.9 


5.6 


io- 


1 ] 


-25.43 


GS 1354-64 


1.0 


3.8 


io- 


1 1 


-64.73 


GRO J0422+32 


0.8 


8.7 


io- 


1 ] 


32.91 


XTE Jill 8+480 


0.7 


1.1 


io- 


10 


48.05 



4. Expected number of microquasar 
events 

In Tab.|2]are given the number of selected neu- 
trino events from each microquasar, applying 
the event selection that optimize the sensitiv- 
ity in Tab. [T] a nd according to the n eutrino 
fluxes given bv lDistefano et al. I d2002l ). A de- 
tailed analysis considering other neutrino pro- 
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Fig. 1. NEMO-km 3 sensitivity to neutrinos from 
microquasars versus source declination, for a live- 
time of 1 year. The worsening of the sensitivity with 
increasing declination is due to the decrease of the 
source visibility. 



duction mod e ls in m icroquasars is reported in 
lAielloetal.l (l2007h . The results in Tab. pre- 
fer to an integration time Af equal to the du- 
ration of the considered burst for the transient 
sources and to 1 year for the steady sources. In 
the same table is given the background (atmo- 
spheric neutrinos + muons) in 1 year of data 
taking. In this analysis, it is assumed that tran- 
sient sources cause one burst per year, i.e. the 
number of source events produced in the inter- 
val Af is relative to 1 year observation time. 

In order to estimate the event rates for 
non-persistent sources it is crucial to know 
their duty cycle. Some of these sources have 
a periodic bursting activi ty: Circinus X-l ha s 
a period of 16.59 days iPreston et al. I d2007l) . 
This means therefore that we expect about 1.5 
events per year. Other transient sources show 
a stochastic bursting activity. For such cases 
it is difficult to give an estimate of the ex- 
pected event rate. For example, during 1994 
GRO J 1 655- 40 had three radio flares, eac h last- 
ing 6 days dHje llming & Rupen 1 1 19951) : dur- 
ing the same year GRS 1915+105 emitted 4 
bursts dRodriguez & Mirabel] 1 19991) . A recent 
study (Ni poti et al. 1 12005) has shown that mi- 
croquasars GRS 1915+105, Cygnus X-3 and 
Scorpius X-l are in flaring mode 21, 10 and 3 
percent of the time, respectively. The possibil- 
ity to integrate over more then one burst could 



therefore help to detect neutrinos from micro- 
quasars. 

The search for neutrino events in coinci- 
dence with microquasar radio outbursts could 
be a tool to reject atmospheric background, re- 
stricting the analysis period to the flare du- 
ration Af. Such an ana lysis technique, al - 
ready used by AMANDA lAckerma nn (2005), 
can improve the detector sensitivity to neutri- 
nos from transient sources. Referring to the 
bursts considered in Tab.|2]and integrating over 
the time interval Af of the bursts, we expect 
an average background of about 10~ 3 events 
(muons) per burst. Summing on all the bursting 
sources, in Tab.|2]we count ~ 0.04 background 
events, which requires about 5 source events 
for a 5cr level detect ion with a 70% probability 
dAhrens et al. | |2004). Tab.|2]shows that we ex- 
pect 3.4 + 9.0 events in the case of a burst from 
each of the bursting microquasars. Therefore, 
a cumulative analysis could provide a possible 
detection of microquasar neutrinos. 

5. Conclusions 

The possibility to detect TeV neutrinos from 
Galactic microquasars with the proposed 
NEMO-km 3 underwater Cerenkov neutrino 
telescope has been investigated. A Monte 
Carlo was carried out to simulate the expected 
neutrino-induced muon fluxes produced by mi- 
croquasars and by atmospheric neutrinos. The 
expected atmospheric muon background was 
also simulated. We computed the detector sen- 
sitivity for each microquasar, optimizing the 
event selection in order to reject the back- 
ground. Finally, we applied the event selec- 
tion and calculated the number of surviving 
events. Our results show that, assuming reason- 
able scenarios for TeV neutrino production, the 
proposed NEMO telescope could identify mi- 
croquasars in a few years of data taking, with a 
discovery potential for at least few cases above 
the 5<r level, or strongly constrain the neutrino 
production models and the source parameters. 
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Source name 


At (days) 


f' h (erg/cm 2 s) 




b 


Steady Sources 


LS 5039 


365 


1.69-10-'" 


0.1 


0.1 


Scorpius X- 1 


365 


6.48- 10 -12 


0.2 


0.1 


SS433 


365 


1.72-10" 9 


76.0 


0.1 


GX 339-4 


365 


1.26-10" 9 


68.0 


0.1 


Cygnus X- 1 


365 


1.88-10' 11 


0.5 


0.1 


Bursting Sources 


XTE J 1748-288 


20 


3.0710 10 


0.8 


0.3 


Cygnus X-3 


3 


4.02- 10~ 9 


0.8 


0.1 


GRO J 165 5 -40 


6 


7.37-10-'° 


0.6 


0.1 


GRS 1915+105 


6 


2.1010 10 


0.1 


<0.1 


Circinus X-l 


4 


1.22-10-'° 


0.1 


0.1 


XTE J 1550-564 


5 


2.00-10-" 


< 0.1 


<0.1 


V4641 Sgr 


0.3 


2.25-10-'° + 3.25 -10- 8 


< 0.1+1.4 


0.1 


GS 1354-64 


2.8 


1. 88-10-" 


< 0.1 


0.1 


GRO J0422+32 


1+20 


2.5110 10 


< 0.1+0.4 


0.1 


XTE Jill 8+480 


30+150 


5.02-10-'° 


1.0+4.8 


0.2 
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